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ABSTRACT 


The  aerobang  maneuver,  one  of  three  types  of  aero-assisted 
orbital  change  methods,  holds  the  possibility  of  reducing  fuel 
consumption  for  orbital  craft  capable  of  atmospheric  entry. 
It  has  been  previously  shown  that  different  types  of  vehicles 
provide  varying  results  over  a  constant  heating  rate 
trajectory.  Further  investigation  into  the  optimization  of 
the  aerobang  maneuver  in  this  thesis  includes  the  effects  of 
using  thrust  vector  control,  the  examination  of  the  effects  of 
increasing  fuel  mass  fraction  to  increase  orbital  inclination 
changes,  and  the  effects  of  that  increase  on  both  angle  of 
attack  and  heating  rate.  The  aerobang  maneuver  is  shown  to  be 
capable  of  significant  changes  in  orbital  inclination  in 
either  a  fixed  heating  rate  or  a  fixed  angle  of  attack  mode 
for  the  Maneuverable  Reentry  Research  Vehicle. 
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I .  INTRODUCTION 


A.  BACKGROUND 


Almost  all  orbiting  spacecraft  have  the  ability  to  correct 
and  change  their  orbits,  either  within  a  given  orbital  plane 
or  by  transferring  to  other  planes  for  better  positioning 
depending  on  the  purpose  of  the  spacecraft.  These  maneuvers 
are  performed  using  the  spacecraft’s  propulsion  and  attitude 
control  systems  and  are  limited  by  the  amount  of  fuel  the 
spacecraft  can  carry  onboard. 

Executing  a  non-coplanar  transfer  reguires  more  fuel  than 
an  orbital  change  within  the  original  plane,  such  as 
increasing  or  decreasing  orbital  velocity  to  change  orbital 
radius  or  altitude.  In  addition  to  velocity  changes,  a  non- 
coplanar  transfer  requires  a  change  in  the  direction  of  travel 
of  the  spacecraft.  Missions  requiring  numerous  plane  changes 
are  greatly  limited  by  the  fuel  capacity  of  the  spacecraft. 

Since  the  early  1960's,  the  possibility  of  using 
aerodynamic  forces  as  a  means  of  controlling  both  velocity  and 
direction  of  a  spacecraft  has  been  studied  [Ref.  1].  With  a 
properly  designed  spacecraft,  the  Earth's  atmosphere  (or  that 
of  any  other  planet)  could  be  used  to  assist  the  propulsive 
force  of  the  spacecraft.  Such  a  spacecraft  would  need  to  be 
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able  to  use  the  lift  forces  generated  to  overcome  the  drag 
present  to  maintain  flight  and  be  constructed  to  withstand  the 
surface  heating  effects  that  would  occur  within  the  atmosphere 
at  near  orbital  speeds.  Figure  1  shows  a  simple  schematic 
drawing  of  what  an  aerodynamically  assisted  orbital  transfer 
would  entail. 


EXOATM.  / 
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Figure  1  Aerodynaalcally  Assisted  Orbital  Transfer 


B.  TYPES  OP  STHER6ETIC  ORBITAL  MAMEUVERS 

Synergetic  maneuvers  are  defined  as  maneuvers  that  use 


both  atmospheric  (aerodynamic)  and  propulsive  forces.  For 


orbital  transfers  between  orbital  planes,  they  have  been 
divided  into  three  categories;  aeroglide,  aerocruise,  and 
aerobang. 

1.  Aeroglide 

As  the  name  implies,  the  aeroglide  maneuver  utilizes 
the  aerodynamic  forces  present  to  create  a  gliding,  unpowered 
trajectory.  While  fuel  is  required  to  deorbit  prior  to  and 
reorbit  after  the  maneuver,  the  aeroglide  itself  relies  only 
on  the  lift  generated  by  the  spacecraft's  interaction  with  the 
atmosphere.  As  a  consequence  of  this,  the  glide  must  not  only 
be  performed  at  a  sufficiently  low  altitude  to  take  advantage 
of  the  more  dense  atmosphere,  but  also  at  velocities  high 
enough  to  maintain  flight.  This  results  in  extremely  high 
surface  heating  rates,  perhaps  beyond  the  capacity  of  the 
spacecraft's  structure  to  absorb  or  dissipate  [Ref.  2].  This 
disadvantage  makes  the  aeroglide  maneuver  less  appealing  for 
non-coplanar  orbital  transfers  by  itself,  although  it  could 
conceivably  be  used  in  conjunction  with  another  technique. 

2 .  Aerocruise 

The  aerocruise  maneuver  uses  the  spacecraft's 
propulsive  force  as  well  as  the  aerodynamic  forces  present  in 
a  glide.  The  drag  encountered  by  the  spacecraft  is  exactly 
compensated  for  by  the  thrust  generated  by  the  spacecraft's 
engines.  This  allows  the  spacecraft  to  maintain  a  constant 
heating  rate  and  avoid  the  overheating  problems  present  with 


the  aerogiide.  This  requires  both  constant  speed  and 
altitude.  The  engine  is  throttled  to  maintain  thrust 
sufficient  to  counteract  drag,  and  the  bank,  angle  is  adjusted 
to  have  a  portion  of  the  lift  oppose  the  cenrrifugal  force, 
thus  maintaining  a  constant  orbital  radius  or  altitude. 
Changes  in  orbital  inclination  come  as  a  result  of  the  lift 
generated  and  will  depend  on  the  lift  to  drag  ratio  (L/D)  and 
the  angle  of  attack  (a)  of  the  vehicle.  Altering  these 
parameters  will  produce  varying  amounts  of  inclination  change 
for  a  given  amount  of  fuel  burned. 

3 .  Aerobang 

This  maneuver  also  uses  a  constant  heating  rate  as  a 
requirement,  allowing  the  altitude  and  velocity  to  vary  during 
the  flight  as  dictated  by  the  control  law  for  the  maneuver, 
discussed  later.  This  allows  a  constant  thrust  to  be  used  and 
allows  the  angle  of  attack  to  vary.  The  thrust  is  set  at  the 
maximum,  reducing  the  time  the  spacecraft  is  in  the  atmosphere 
compared  to  the  aerocruise,  thus  reducing  the  integrated 
heating  effects  on  the  spacecraft.  It  has  previously  been 
shown  that,  for  certain  velocities,  altitudes,  and  heating 
rates,  an  aerobang  maneuver  will  yield  a  greater  orbital 
inclination  change  than  an  aerocruise  maneuver,  for  the  same 
fuel  expenditure  [Ref.  3]. 


C.  SCOPE  OF  THESIS 


This  thesis  will  further  examine  the  aerobang  non-copianar 
orbital  transfer  maneuver.  The  initial  results  showed  that 
the  maneuver  generated  smaller  than  expected  inclination 


changes  for 

the 

given  two  percent 

mass 

fraction 

of 

fuel 

consumed. 

The 

thrust  acted  in 

the 

direction 

of 

the 

spacecraft ' s 

body  and  thus  only  a  portion 

(depending 

upon 

the 

angle  of  attack)  acted  to  help  the  lift,  which  produced  the 
inclination  change.  The  effects  of  thrust  vector  control 
(TVC)  are  included  here  to  allow  the  thrust  to  contribute 
further  to  the  lift,  and  thus  increase  the  inclination  change 
while  maintaining  an  angle  of  attack  within  the  vehicle's 
capability. 

The  aerobang  maneuver  is  limited  by  the  amount  of  fuel 
allotted  for  each  atmospheric  pass.  Obviously,  the  greater 
the  fuel  used,  the  longer  the  maneuver  can  be  performed, 
subject  to  heating  and  altitude  restraints.  The  initial 
results  [Ref.  3]  were  based  upon  a  0.02  mass  fraction  because 
the  aerobang  maneuver  was  unable  to  maintain  a  constant 
heating  rate  profile  beyond  that  point  for  some  conditions  of 
altitude  and  velocity.  The  inclination  changes  generated  with 
this  restriction  (less  than  one  degree  for  the  most  part)  is 
not  useful  for  many  orbiting  spacecraft  applications.  By 
expanding  the  allowable  mass  fraction  of  the  fuel  consumed  for 
the  maneuver  to  values  up  to  twenty  percent,  larger,  more 
useful  inclination  changes  should  be  able  to  be  produced. 
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In  addition  to  the  two  limiting  heating  rates  examined  in 
Reference  3,  the  effects  of  both  aerocruise  and  aerobang  over 
the  entire  range  of  heating  rates  are  studied.  The  intent  is 
to  determine  the  shape  of  the  heating  rate  profile  for  both 
maneuvers  and  determine  more  precisely  when  aerobang  is  a 
superior  method  of  orbital  transfer.  This  study  will  also 
observe  the  behavior  of  the  vehicle's  angle  of  attack,  during 
extended  constant  heating  rate  maneuvers  as  well  as  fixing  the 
angle  of  attack  for  other  maneuvers  and  observing  its  effect 
on  the  heating  rate. 

D.  FLIGHT  VEHICLE 

The  vehicle  used  for  this  study  is  the  Maneuverable 
Reentry  Research  Vehicle  (MRRV).  The  MRRV  was  designed  as  a 
hypersonic  research  vehicle  to  be  launched  by  the  space 
shuttle  and  capable  of  a  powered  re-entry  and  flight  through 
the  atmosphere  to  demonstrate  the  feasibility  of  aerocruise 
maneuvers  [Ref.  4].  The  MRRV  has  also  been  used  as  a  test 
vehicle  for  aerobang  simulations  [Ref.  3].  Because  of  this, 
and  for  consistency,  it  is  used  here  as  well. 

As  depicted  in  Figure  2  [Ref.  3],  the  MRRV  is  a  winged 
body  approximately  7.6  meters  long  with  a  four  meter  wingspan. 
The  effective  surface  area  is  11.61  with  an  initial  mass  of 
4898  kg,  of  which  2588  kg  is  fuel.  The  MRRV  is  powered  by 


three  Marquardt  R-40-B  rocket  motors  that  provide  a  total 
thrust  of  14679  Newtons  with  a  specific  impulse  (I^p)  of  295 
seconds . 


Figure  2  The  MRRV 


Figure  3  shows  that  the  MRRV  has  a  nominal  lift  to  drag 
ratio  (L/D)  of  2.3  at  its  optimal  angle  of  attack  of  14 
degrees.  This  is  based  upon  wind  tunnel  data  gathered  out  to 
an  angle  of  attack  of  between  30  and  40  degrees  [Ref.  4].  No 
data  was  available  for  higher  values,  implying  that  this  is 
the  upper  limit  for  angle  of  attack  for  this  vehicle. 


Figure  3  C^/C,  and  Angle  of  Attack  for 

the  NRRV 


E.  APPLICATIONS 

Initially,  it  would  appear  that  the  results  described  in 
^  (and  modified  later  in  this  thesis  for  the  2% 
case)  may  not  be  useful  for  most  orbital  transfer  applications 
since  they  amount  to  less  than  a  one  degree  change  in 
inclination,  occurring  over  a  period  of  under  20  seconds. 
There  are,  however,  possible  applications  and  factors  to  take 
into  account  when  studying  these  results. 

An  inclination  change  of  even  a  degree  or  so  can 
correspond  to  a  difference  of  over  a  hundred  kilometers  when 
the  angle  is  projected  up  to  the  orbital  altitudes  used  here 
and  higher  as  well.  For  example,  for  a  one  degree  change  in 
inclination  at  a  radius  of  6447  km,  the  spacecraft  will  move 
112.5  km.  In  a  military  application  involving  pursued  and 
pursuing  vehicles  in  orbit,  such  as  anti-satellite  warfare,  a 
quick  aerobang  maneuver  resulting  in  even  a  small  amount  of 


■9) 

♦ 


8 


•  •••••• 


•  • 


orbital  transfer  may  be  capable  of  preventing  a  vehicle  from 
being  detected  or  destroyed. 

Even  for  civilian  applications,  a  small  aerobang  maneuver 
may  be  preferred  to  a  longer  one  to  minimize  the  heat  load 
imparted  to  a  vehicle  that  may  not  have  the  thermal  protection 
system  to  adeguately  protect  against  a  longer  atmospheric 
flight  that  would  yield  a  greater  inclination  change.  By 
making  multiple,  low  fuel  passes  and  reorbiting  between 
orbital  changes  to  radiate  the  heat  to  space,  the  same  total 
inclination  change  can  be  made  if  time  is  not  of  the  essence. 


II.  EQUATIONS  OF  MOTION 


♦ 

The  motion  of  the  aerodynamic  spacecraft  within  its  orbit 
is  defined  by  six  variables  for  both  position  and  velocity.  • 

Position  is  defined  by  the  standard  spherical  coordinates;  r 
(radius),  0  (right  ascension),  and  <p  (declination).  The  other 
three  variables  are  v  (velocity),  y  (flight  path  angle),  and  * 

♦  (heading  angle).  These  variables  are  graphically  shown  in 
Figure  4  below. 


The  flight  path  angle  is  described  as  the  angle  of  the 
velocity  vector  with  relation  to  the  local  horizontal  plane 
tangent  to  the  orbit  and  the  heading  angle  is  the  direction 
angle  within  that  plane  relative  to  the  local  inertial 
latitude . These  equations  were  derived  [Ref.  3)  assuming  a 
spherical,  non-rotating  Earth .  Atmospheric  motion  (winds)  and 
the  oblateness  of  the  Earth  were  also  neglected.  The  method 
of  derivation  and  the  resulting  equations  that  were  used  in 
this  study  are  summarized  in  this  chapter. 

A.  SPACECRAFT  POSITION 

The  derivation  of  the  equations  of  motion  requires  the  use 
of  three  coordinate  systems.  As  shown  in  Figure  4,  the 
position  of  the  vehicle  can  be  expressed  in  the  inertial  XYZ 
system,  the  orbital  xyz  system,  also  centered  at  the  Earth's 
center,  with  the  positive  x  axis  pointed  to  the  orbiting 
spacecraft,  as  well  as  the  x'y'z'  system  which  is  centered  on 
the  spacecraft  with  the  positive  y'  direction  in  line  with  the 
velocity  vector  and  the  x'  axis  within  the  xy  plane. 

Direction  cosine  matrices  (DCM)  are  required  to  transform 
the  various  coordinate  systems  into  one  to  be  used  for  the 
expression  of  the  spacecraft's  position  and  velocity.  In 
order  to  get  from  the  inertial  XYZ  system  to  the  orbital  xyz 
system,  a  body  3-2  rotation  through  angles  0  and  (p  is 
required.  To  obtain  the  xyz  coordinates  from  the  x'y'z' 
system  a  3-1  rotation  is  used  through  angles  y  and 
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By  equating  the  velocity  vector  to  the  time  derivative  of 
the  position  vector  for  the  spacecraft,  the  kinematics  can  be 
expressed  as  follows: 


^  =  Vsiny 
dt 

(2.1) 

dd  _  '/cosycostlf 
dt  rcos^ 

(2.2) 

d<t>  _  l/cosYsirn|f 

(2.3) 

dt  r 

The 

other  three  equations  are  based  on  the  balance  of 

forces 

on  the  spacecraft,  and  their 

relation  to  the 

accelerations  and  masses  (F*ma).  The  forces 

are  first  derived 

in  the  spacecraft's  frame  of  reference,  the 

x'y'z'  system,  as 

shown  in 

Figure  5. 

The 

forces  in  the  x',  y'  and  z'  directions  are: 

[  Tsin(a*€)  *L]  cosa 

(2.4> 

F^i=Tcos{a*€)  -D 

(2.5) 

F^,=  [Tsin{a*€)  ]  -D 

(2.6) 

Again,  the  conversion  to  the  rotating  reference,  the  xyz 
system,  is  done  using  the  appropriate  DCM.  The  resulting 
expressions  for  the  forces  on  the  spacecraft  are  then  set 


equal  to  the  product  of  the  mass  and  the  accelerations 
(calculated  by  taking  the  derivative  of  velocity  vector)  in 
each  of  the  three  direction  components  (x,  y,  and  2)  and  are 
reduced  to  the  following: 

[rcos  (o*€) -D]  -/njgsinY=/77^^  {2.1) 

dy  \Tsin(a*(.)  *L]coso-m^gcosy  o  > 

- - - COSY  (2.8) 

dt  r 

dt  m^cosy  r 

The  components  of  the  spacecraft's  acceleration  in  the 
tangential,  normal,  and  binormal  directions  can  be  expressed 
independently.  This  will  allow  the  final  t..ree  equations  of 
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motion  to  be  expressed  in  a  more  compact  form.  For  a 
perfectly  circular  orbit,  the  normal  (to  the  tangential 
velocity  and  acceleration  vectors)  direction  is  also  the 
radial  direction,  pointing  to  the  center  of  the  Earth.  Hence, 
the  notation  for  this  component. 

A,=  [7-sin(a^€)^f.:coso  ^2.10) 


Tcos(a^€>-D  ^2.11) 

[rsin(a>€,^f.;sino 

The  final  three  eguations  of  motion  then  become: 

^=A^-gsiny  (2.13) 

V-^  =A,+ ( — -g)  cosy  (2.14) 

at  r 

+-^cosYCOSi|»tan<|>  (2.15) 

dc  cosy  r 


T" 


;4 


j 


•  ( 


B.  ORBITAL  INCLINATION 

The  measure  of  the  effectiveness  of  the  orbital  plane 
change  maneuver  is  the  change  in  the  orbit's  inclination. 
Since  the  equations  of  motion  use  the  spherical  coordinates  to 
describe  the  spacecraft's  motion,  the  inclination  must  also  be 
expressed  in  these  coordinates.  As  shown  in  Reference  3,  the 
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angular  momentum  vector  will  be  perpendicular  to  the  orbital 
plane  so  that  the  angle  between  that  vector  and  the  unit 
vector,  k,  defining  the  positive  Z  direction  will  be  the 
orbit's  inclination.  The  magnitude  of  the  angular  momentum 
vector  is  shown  to  be  Vrcosy  while  the  Z  component  of  the 
angular  momentum  is  VrcosYcosiji .  Since 

h-k=\f\cosi  (2.16) 

then 


l''rcosYCOstltcos<J>=|/ilcos  i 
and  so  it  follows  that 

cosi  =cos»|;cos4> 


(2.17) 


(2.18) 


giving  a  means  of  simply  calculating  the  orbital  inclination 
for  any  plane  change  maneuver. 


III.  HEATING  RATE  AMD  CONTROL  LANS 


The  control  parameters  for  the  aerobang  maneuver  are  the 
bank  angle,  the  angle  of  attack,  and  the  thrust.  The  angle  of 
bank  is  set  to  90  degrees,  the  thrust  is  set  to  the  maximum, 
and  the  angle  of  attack  is  adjusted  to  maintain  a  constant 
heating  rate  on  the  spacecraft  during  the  maneuver.  In  order 
to  implement  a  control  law  for  the  aerobang  orbital  transfer 
maneuver,  the  modeling  of  the  atmosphere  and  resulting  heating 
rate  must  first  be  developed. 

A.  ATMOSPHERIC  MODEL 

The  first  step  in  the  modeling  of  the  heating  rate 
behavior  is  to  model  the  atmosphere  in  the  region  in  which  the 
maneuvers  will  take  place.  For  the  aerobang  maneuvers 
examined  in  this  and  previous  studies,  the  region  between  the 
altitudes  of  50  and  120  km  was  used.  The  density  of  the  air 
at  these  altitudes  was  modeled  using  the  U.S.  Standard 
Atmosphere,  1976,  [Ref.  5,  p.  7)  and  is  approximated  by  a 
local  exponential  model 


where 


p,  =  3.0968x10'*  kg/m^ 


B.  HEATIMG  RATE  MODEL 


R  =  1.41xl0'‘  m*'' 
=  6438  km 


Atmospheric  density  (and  therefore  altitude,  as  shown 
above)  along  with  the  vehicle's  velocity  combine  to  give  a 
basic  relationship  for  the  heating  rate.  As  shown  in 
Reference  5,  that  relationship  is: 

Q^kp-''V'"  (3.2) 

with  k,  n,  and  m  as  constants.  This  eguation,  combined  with 
the  equations  of  motion  can  be  used  to  derive  the  angle  of 
attack  required  to  fly  at  a  constant  heating  rate  during  the 
aerobang  maneuver . 

By  combining  Equations  3.1  and  3.2,  the  relationship 
between  density,  heating  rate,  and  velocity  becomes  more 
apparent . 


r 

v=\- 


P  (r-r,;  I 


(3.3) 


,  J 

In  order  to  show  the  effect  of  the  heating  rate  model  on  the 
equations  of  motion,  the  derivative  of  Equation  3.3  is  taken, 
and  simplified  to  give: 


dt  i  m  ]  dt 


(3.4) 


Using  Equation  2.1,  the  expression  for  dr/dt  can  be  used  to 
express  dV/dt  as  a  function  of  V  and  y ,  the  flight  path  angle. 
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V^s iny 


(3.5) 


dVJ 
dt  m 


Using  Equations  2.11  and  2.13  gives  for  the  equation  along  the 
path, 


dV _  Tcos  (g ■»•€)  -D 
dt 


-gsiny 


(3.6) 


C.  DETERMINING  THE  ANGLE  OF  ATTACK 

Equating  Equations  3.5  and  3.6  gives  an  expression  for 
calculating  the  angle  of  attack  required  for  the  aerobang 
maneuver : 

0  =  7’cos  (o^c)  -D-m  siny  (g*-^  V^)  (3.7) 

in 

The  drag  term,  D,  is  also  an  function  of  a  so  the  solution  for 
alpha  is  given  by  the  transcendental  equation  above. 

In  order  to  solve  for  an  angle  of  attack  to  be  used 
throughout  the  maneuver,  an  initial  value  must  be  selected  and 
the  above  equation  (3.7)  must  be  iterated  to  solve  for  a  at 
each  successive  time  interval.  As  shown  in  Reference  3, 
Newton's  Method  of  convergence  was  used  to  arrive  at  an  angle 
of  attack  for  this  situation.  Initially  guessing  a  involved 
using  a  binomial  curve  fit  for  the  drag  versus  angle  of  attack 
data.  Assuming  minimal  values  for  the  tangential  acceleration 
and  angle  of  attack  (so  that  the  cosine  is  nearly  equal  to 
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one)  and  no  thrust  vector  angle  in  Equation  2.11,  then  the 
thrust  is  equal  to  the  drag  (or  the  binomial  expansion  curve 
fit  for  drag) , 


r 


(3.8) 


where  C^,  Cj,  and  are  the  coefficients  for  the  zeroth,  first 
and  second  order  terms,  respectively,  in  the  binomial  curve 
fit  expression  [Ref.  3:  p.  13].  The  quadratic  formula  gives 
a  method  for  calculating  the  first  guess  for  a. 


a  = 


1 

-cl 

1  * 

N 

2C, 


(3.9) 


This  technique  is  used  in  the  computer  program  to  simulate 
the  aerobang  maneuver  to  get  a  satisfactory  approximation  for 
the  varying  angle  of  attack. 


D.  THE  HEATING  RATE  PROFILE 

The  previous  thesis  examined  the  aerocruise  and  aerobang 
at  only  two  heating  rate  values.  It  was  assumed  that  the 
resulting  regions  of  superior  performance  for  both  subcircular 
and  supercircular  speeds  held  for  all  heating  rates  in  between 
[Ref.  3:  p.  63]. 

Prior  to  a  closer  examination  of  the  behavior  of  flight 
parameters  for  the  aerobang  maneuver,  this  assumption  was 
validated  by  determining  the  velocities  at  which  aerocruise 


lost  its  superiority  to  aerobang  in  terms  of  the  inclination 


changes  produced  during  the  2%  mass  fraction  runs.  These 
results,  shown  in  Figure  6,  do  show  that  at  very  near  circular 
speeds,  the  aerocruise  maneuver  did  produce  better  results, 
but  that  the  profile  is  not  entirely  symmetrical  about  the 
circular  velocity.  In  general  it  was  true  that,  at  lower 
heating  rates,  the  aerocruise  method  was  better  for  a  much 
smaller  range  of  velocities. 
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Figure  6  Aerobang  and  Aerocruise  Efficiency 
Velocities  and  Beating  Rates 
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IV.  SIMULATION  PROGRAM 


In  order  to  evaluate  the  potential  of  the  aerobang  orbital 
plane  change  maneuver,  a  FORTRAN  program,  previously  developed 
[Ref.  3],  was  modified  and  used  to  simulate  the  action  of  the 
spacecraft  while  in  the  atmosphere.  The  modifications  to  the 
original  program  included  the  addition  of  £,  the  thrust  vector 
angle,  the  addition  of  a  constant  angle  of  attack,  option,  and 
the  elimination  of  the  aerocruise  option  and  the  associated 
calculations . 

The  program,  listed  in  Appendix  A,  was  written  in  a 
modular  format,  using  several  subroutines  for  calculations, 
integration  of  the  equations  of  motion,  and  output.  The  main 
program  calls  each  one  and  controls  their  use  in  the  program. 
For  the  heating  rate  versus  velocity  simulations  used  in 
Chapter  III,  the  original  program  [Ref  3;  App.  A]  was  used  to 
produce  the  aerocruise  data  shown  in  Figure  6. 

A.  THE  MAIN  PROGRAM 

The  main  portion  of  the  program  links  the  subroutines 
together  and  controls  the  iteration  of  the  solution.  The  two 
input  files  and  one  output  file  (also  in  Appendix  A)  are 
opened  first.  The  use  of  the  two  input  files  (AERO.DAT  and 
DATA.DAT)  minimized  the  changes  that  had  to  be  made  to  the 
program  and  subroutines  themselves.  The  AERO.DAT  file 
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contains  the  information  about  the  spacecraft  and  the  DATA.DAT 
file  contains  the  information  concerning  the  flight  profile 
and  the  printing  interval.  The  main  program  also  calculates 
the  mass  changes  throughout  the  simulation  and  stops  the 
program  when  the  final  mass  is  reached  or  if  the  time  extends 
beyond  the  maximum  allowed.  A  flow  chart  for  the  main  program 
appears  in  Figure  7  [Ref.  3]. 

B .  THE  SUBROUTINES 

1.  The  CHTRL  Subroutine 

The  control  subroutine  is  the  major  subroutine, 
calculating  the  atmospheric  and  heating  rate  models,  as  well 
as  the  angle  of  attack  for  the  aerobang  control  itself.  The 
angle  of  attack  is  computed  using  the  Newton  approximation 
routine  as  discussed  in  Chapter  II  to  simulate  the  modulation 
of  the  angle  of  attack  for  the  aerobang  maneuver. 
Alternatively,  the  angle  of  attack  can  be  fixed  and  allowed  to 
change  the  heating  rate  accordingly.  In  either  case,  the  lift 
and  drag  forces  are  calculated  here  for  use  later  in  the 
equations  of  motion.  Figure  8  shows  the  flow  chart  for  this 
subroutine . 

2.  The  ACEL  Subroutine 

This  subroutine  computes  the  three  accelerations, 
defined  by  Equations  2,10,  2.11,  and  2.12.  These  values  are 
used  in  the  next  subroutine  to  evaluate  the  motion  of  the 
spacecraft . 
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Figure  7  Flow  Chart  for  the  Main  Program 
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3.  The  ORB  Subroutine 


In  this  portion  of  the  program,  the  six  equations  of 
motion  (Equations  2.1,  2.2,  2.3,  2.13,  2.14,  and  2.15)  are 
defined  in  terms  of  the  six  variables  defining  the  motion  of 
the  spacecraft;  given  as  a  vector,  X(6),  along  with  their 
derivatives,  XD0T(6).  These  variables  are  used  in  the 
subroutine  to  compute  the  motion  of  the  spacecraft  at  each 
time  iteration.  The  equations  are  integrated  in  a  Runge-Kutta 
routine,  described  later. 

4.  The  WRT  and  HDR  Subroutines 

These  subroutines  control  the  output  of  the  program. 
The  HDR  (header)  routine  creates  a  header  in  the  output  file 
that  gives  the  initial  conditions  of  the  program  and  the 
headers  for  the  columns  of  data  that  the  program  computes. 
The  WRT  (write)  subroutine  is  called  from  the  main  program  and 
used  on  each  time  iteration  to  write  the  calculated  data  into 
the  output  file.  It  is  also  called  at  the  beginning  and  end 
of  the  iterations  to  put  both  the  initial  and  final  data  rows 
into  the  file. 

5.  The  RK4  Subroutine 

This  subroutine  is  the  fourth  order  Runge  Kutta 
routine  used  to  integrate  the  equations  of  motion.  It  was 
written  by  Professor  I.  M.  Ross  and  is  a  part  of  the  public 
collection  of  subroutines  at  the  Naval  Postgraduate  School. 
Here,  the  time  iterations  are  counted  and  used  to  create  the 


multiple  integration  computations  and  send  the  results  to  the 
main  program  for  writing  into  the  output  file. 

C.  PROGRAM  VERIFICATION 

Since  this  program  is  based  on  an  earlier  version,  the 
validation  of  ;his  program  involved  using  it  to  verify  results 
obtained  in  the  previous  study.  By  setting  the  thrust  vector 
angle  egual  to  zero  in  the  AERO.DAT  file's  last  entry,  the 
exact  same  results  should  be  obtained  for  similar  initial 
conditions  used  previously. 

The  comparison  of  the  aerobang  and  aerocruise  maneuvers  to 
verify  the  applicability  of  aerobang  over  the  original  heating 
rate  profile  provided  an  opportunity  to  verify  the  program's 
proper  operation.  By  using  both  programs  to  calculate  results 
for  the  aerobang  cases  and  obtaining  the  same  results,  the  new 
program  could  then  be  used  on  its  own  to  incorporate  the 
effects  of  thrust  vector  control  and  serve  as  a  tool  for 
further  analysis  of  the  aerobang  maneuver. 


V.  AEROBANG  PERFORMANCE  ANALYSIS 


As  shown  in  Chapter  III  and  previous  studies  [Ref.  3], 
there  exist  certain  conditions  in  which  the  aerobang  orbital 
plane  change  is  the  most  efficient  as  determined  by  the  amount 
of  orbital  inclination  change  that  can  be  achieved  for  a  given 
amount  of  fuel  expended.  Having  determined  the  conditions 
necessary  for  successful  aerobang  maneuvers,  the  next  step  is 
to  determine  the  ways  in  which  the  aerobang  maneuver  can  be 
improved  and  best  utilized  to  achieve  a  given  mission 
criterion.  The  effects  of  adding  a  thrust  vector  control 
(TVC)  system  to  improve  upon  previously  obtained  optimization 
results  was  examined  first.  Further  simulations  focused  on 
the  effects  of  increasing  the  mass  fraction  of  fuel  burned  to 
allow  for  more  useful  inclination  changes  (on  the  order  of  ten 
or  more  degrees  per  maneuver)  and  the  advantages  and 
disadvantages  of  fixing  the  angle  of  attack  rather  than  the 
heating  rate  for  a  specific  maneuver. 

In  order  to  remain  consistent  with  the  previous  study 
[Ref.  3]  and  to  be  able  to  compare  the  previous  results  with 
those  obtained  and  given  here,  a  review  of  those  results  is 
necessary . 

The  previous  thesis  looked  at  both  subcircular  and 
supercircular  velocities  for  both  aerocruise  and  aerobang 
maneuvers  for  two  distinct  heating  rates.  The  effects  of  the 
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heating  rate  values  over  the  entire  range  rather  than  just  the 
two  endpoints  on  the  viability  of  both  subcircular  and 
supercircular  results  were  discussed  in  Chapter  III.  These 
results  were  based  on  a  mass  fraction  of  fuel  expended  of  only 
two  percent. 

For  the  higher  heating  rate  (1.42x10'^  W/m^),  the  best 
aerobang  performance  results,  in  terms  of  inclination  changes, 
occurred  at  a  velocity  factor  (k=V^/V)  of  0.98  for  the 
supercircular  case  and  1.02  for  the  subcircular  case,  with 
inclination  changes  of  0.554°  and  0.789°,  respectively. 

A.  THRUST  VECTOR  CONTROL 

1 .  Subcircular  Case 

In  order  to  determine  the  overall  effects  of  thrust 
vector  control  on  the  aerobang  maneuver,  the  computer 
simulations  were  run  for  a  range  of  c  from  -80  to  80  degrees. 
While  these  angles  may  not  presently  be  feasible  from  a 
structural  or  control  point  of  view,  they  were  included  for 
completeness.  The  simulations  were  run  at  an  initial  radial 
distance  of  6447  km,  corresponding  to  an  altitude  of  77  km, 
and  a  heating  rate  of  1.42  MW/m^.  As  shown  in  Figure  9  below, 
the  most  efficient  value  for  e  fell  between  10°  and  16°,  giving 
an  inclination  change  of  0.8°  over  that  range  of  vectoring 
angles.  For  negative  values  of  thrust  vector  angles,  the 
inclination  changes  were  not  as  large  as  those  at  zero  and 
positive  angles  since  negative  values  of  e  will  have  a 
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component  in  the  direction  directly  opposing  the  lift  which  is 
responsible  for  the  change  in  inclination. 


Figure  9  Thrust  Vector  Angle  Effects  -  Subcircular  Case 
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2 .  Supercircular  Case 


For  k  =  0.98,  at  the  same  altitude  and  heating  rate, 
similar  results  were  obtained  for  the  addition  of  thrust 
vectoring  for  this  case.  As  shown  in  Figure  10,  the  best 
inclination  change  occurred  with  a  thrust  vector  angle  of  15 
to  20  degrees,  resulting  in  an  inclination  change  of  0.570°. 

At  higher  c,  the  angles  of  attack  decreased  as  the  spacecraft 
moved  faster  to  higher  less  dense  air  with  a  subsequent  loss 
of  lift. 
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Figure  10  Thrust  Vector  Angle  Effects  - 
Supercircular  Case 
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The  shape  of  the  curve  over  the  broad  range  of  thrust  vector 
angles  is  again  indicative  of  the  effects  of  extreme  values  of 
£  on  flight  performance. 


B.  INCREASING  FUEL  MASS  FRACTION 

It  would  be  logical  to  conclude  that  by  increasing  the 
amount  of  fuel  consumed  for  an  aeiobang  maneuver,  a  greater 
change  in  inclination  would  be  achieved.  By  establishing  a 
goal  of  achieving  an  inclination  change  of  10°  or  more,  the 
effects  of  increasing  the  allowable  mass  fraction  were 
investigated.  These  simulations  were  run  for  both  subcircular 
and  supercircular  orbits  (k=1.02  and  0.98),  at  an  initial 
altitude  of  77  km  (R=6447  km),  and  restricted  to  the  same 
heating  rate  as  was  used  in  the  previous  2%  fuel  mass  fraction 
runs . 


1 .  Subcircular  Case 

a.  Constant  Heating  Rate 

Initially,  the  increase  of  mass  fraction  was  tried 
for  the  constant  heating  rate  case  usually  prescribed  for  the 
aerobang  maneuver.  As  the  mass  fraction  was  increased  from 
the  original  2%  up  to  20%,  it  became  apparent  that  the 
technique  of  modulating  the  angle  of  attack  to  maintain  a 

m 

constant  heating  rate  profile  was  difficult  to  implement. 

Figure  11  shows  the  expected  effect  of  increasing  mass 
fraction  on  aerobang  maneuvers.  The  inclination  change  is 


very  impressive,  between  25°  and  28°,  depending  upon  the  thrust 
vector  angle.  In  these  simulations,  five  values  of  £  (0,  10, 
20,  30,  and  40°)  were  used  to  show  the  effects  of  TVC  on  these 
extended  runs.  The  thrust  vector  angle  labels  for  this  figure 
refer  to  the  right  hand  side  of  the  curves;  zero  at  the  top, 
40  at  the  bottom.  The  smaller  values  of  e,  gave  the  greater 
inclination  changes.  While  this  seems  to  contradict  the  2% 
mass  fraction  results,  the  data  for  time  values  less  than  60 
seconds  revealed  that  the  behavior  of  the  inclination  followed 
the  behavior  shown  earlier  for  the  19  second,  2%  mass  fraction 
case.  The  e  =  10°  curve  starts  out  above  the  0°  curve,  but  the 
curves  cross  at  56  seconds  and  the  0°  curve  finishes  slightly 
above  the  10°  curve.  The  two  curves  are  nearly 
indistinguishable  from  each  other,  but  the  results  for  the  2% 
case  (the  slight  advantage  of  the  ten  degree  thrust  angle  with 
the  three  other  cases  even  lower)  are  shown  in  Figure  12  and 
are  consistent  with  the  results  for  the  longer  times  of  Figure 
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Figure  11  Inclination  For  Constant  Heating  Rate 


Figure  12  Magnification  of  the  Constant  Heating  Rate 
Inclination 


This  superior  performance  in  orbit  transfer  comes  with  a 
penalty.  Figure  13  shows  the  change  in  angle  of  attack  over 
the  extended  simulations.  Values  of  a  rise  to  between  65  and 
70  degrees,  depending  again  on  the  thrust  vectoring,  and 
decrease  to  about  53  or  54°.  The  zero  thrust  angle  case  had 
the  highest  angles  of  attack  throughout  the  flight.  The 
addition  of  thrust  vectoring  caused  the  angle  of  attack  to 
decrease.  This  reduction  in  a  reduced  the  lift  generated  and 
the  inclination  changes  as  well. 


The  aerodynamic  data  provided  for  the  MRRV,  mentioned  in 
Chapter  I,  only  covered  angles  of  attack  of  up  to  40°.  It  can 
be  assumed,  therefore,  that  the  vehicle  is  probably  not 
capable  of  flying  at  a's  greater  than  40°.  If  that  is  the 
case,  then  a  vehicle  of  the  same  surface  area,  initial  mass, 
fuel  burning  rate,  and  thrust  would  not  meet  these  high  angle 
of  attack  results,  if  they  are  valid,  while  maintaining  the 
constant  heating  rate  profile  simulated  in  this  study. 
b.  Constant  Angle  Of  Attack 

At  this  point,  a  departure  from  the  standard 
aerobang  maneuver  theory  was  attempted  to  try  to  avoid  the 
angle  of  attack  problems  described  above.  Rather  than 
allowing  the  angle  of  attack  to  vary  and  allowing  the  constant 
heating  rate  control  the  maneuver,  the  angle  of  attack  was 
fixed  and  the  heating  rate  was  allowed  to  vary  to  evaluate  the 
possible  inclination  changes  while  observing  the  heating  rate 
to  see  if  it  would  exceed  the  desired  heating  rate  and,  if  so, 
by  how  much. 

An  inclination  change  of  approximately  10°  was  still  a 
goal  and  several  values  of  a  were  tested  to  see  if  this  goal 
could  be  met.  The  graphs  in  Figures  14,  15,  and  16  show  the 
results  of  these  runs.  Angles  of  attack  of  30°,  35°,  and  40°, 
respectively,  were  selected  to  get  the  desired  ten  degree 
inclination  change.  Again,  the  same  five  values  of  thrust 
vectoring  are  included  to  show  their  effects.  As  a  increases. 


the  inclination  changes  rise  from  values  of  6.5°  to  9.5°  at  a 


=  30°  to  values  of  12°  to  17°  at  a  =  40°.  In  this  case,  the  • 

increasing  thrust  angle  with  the  fixed  angle  of  attack 

resulted  in  increasing  inclination  changes  throughout  the 

runs.  With  a  fixed,  its  contribution  to  the  vehicle's  lift  is  • 

fixed  and  the  changes  in  inclination  for  this  case  are  made  by 

the  contribution  of  the  thrust  vectoring  to  the  lift 

generated.  • 


Figure  14  Inclination  For  30°  Angle  of  Attack 
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The  result  of  more  concern  is  that  of  the  heating  rates 
that  existed  as  a  result  of  the  fixed  a  flight  regime. 
Figures  17,  18,  and  19  show  the  heating  rates  obtained  for 
these  simulations.  They  vary  from  1.8  x  lo*  W/m^  on  up  to  2.2 
X  10*  W/m^,  both  significantly  above,  almost  up  to  twice  the 
value  used  for  all  previous  runs.  A  vehicle  used  to  perform 
these  fixed  a  aerobang  maneuvers  would,  therefore,  need  to 
have  a  thermal  protection  system  capable  of  withstanding  these 
heating  rates. 


•J 
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2 .  Supercircular  Case 

a.  Constant  Heating  Rate 

Supercircular  aerobang  maneuver  simulations  using 
20%  fuel  mass  fraction  were  attempted  over  the  same  range  of 
angles  of  attack  and  thrust  vector  angles  as  for  the 
subcircular  case.  Under  otherwise  identical  conditions,  the 
supercircular  simulations  were  unable  to  generate  any  results 
that  were  even  close  to  those  of  the  subcircular  case.  The 
inclination  changes  were  no  greater  than  0.865°  for  an  e  of  10° 
and  the  burn  time  was  no  more  than  37  seconds,  approximately 
4%  mass  fraction  as  opposed  to  the  desired  20%  set  in  the 
simulation  program.  The  combination  of  decreasing  a  and 
increasing  altitude  at  the  supercircular  velocity  caused  the 
vehicle  to  lose  its  lifting  force,  resulting  in  the  shortened 
flight.  Compared  to  the  subcircular  case,  the  supercircular 
aerobang  maneuvers  proved  to  be  ineffective  in  producing 
substantial  inclination  changes  (on  the  order  of  10°)  for  any 
kind  of  significant  increases  in  fuel  mass  fraction,  although 
they  may  still  be  viable  for  short,  quick  maneuvers.  Figures 
20  and  21  show  the  inclination  and  angle  of  attack  results  for 
this  case. 
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b.  Constant  Angle  of  Attack 

This  case  proved  to  be  more  successful  since  the 
angle  of  attack  (and  therefore  the  lift  generated)  never 
decreased  to  the  values  near  zero  as  it  did  in  the  constant 
heating  rate  case.  Each  aerobang  simulation  was  completed  for 
the  entire  20%  burn  and  significant  inclination  changes  were 
observed.  Figures  22,  23,  and  24  show  these  results. 


Figure  22  Inclination  for  30°  Angle  of  Attack 
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As  in  the  subcircular  case,  the  inclinations  here 
increased  with  higher  angle  of  attack  and  with  increased 
thrust  vector  angles.  The  maximum  values,  however,  were  not 
as  large  since  the  supercircular  trajectories  caused  the 
altitude  of  the  vehicle  to  increase  and  the  resulting 
atmospheric  density  to  decrease,  thus  reducing  the  lift 
capability . 

The  heating  rates  for  supercircular  flights,  however, 
showed  an  inherent  advantage  with  this  case.  The  higher 
altitudes  had  a  gieater  effect  on  heating  rates  than  the 
increased  velocities  did,  resulting  in  decreasing  values  for 
heating  rate.  As  shown  in  Figures  25,  26,  and  27;  the  heating 
rates  for  all  cases  started  with  the  same  value  of  1.45  x  lO* 
W/m^.  The  final  values  depend  on  both  angle  of  attack  and 
thrust  vector  angle.  The  higher  angle  of  attack  cases 
resulted  in  less  pronounced  decreases  in  heating  rate  value 
from  start  to  finish.  The  thrust  vector  angles  in  these  cases 
actually  contributed  to  higher  heating  rates  as  well,  due  to 
the  slight  increase  in  lift  and  velocity  that  they  provided. 
In  particular,  as  seen  in  Figure  27,  a  limit  is  reached  in  the 
40°  thrust  vector  case,  where  the  heating  rate  reaches  a 
minimum  limit  and  then  begins  to  increase  again. 


HEATING  RATE  (W/M2  X  E 


Figure  25  Heating  Rates  For  30°  Angle  of  Attack 


Figure  26  Heating  Rates  For  35°  Angle  of  Attack 


VI.  CONCLUSIOMS  AMD  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  aerobang  orbital  plane  change  maneuver  has  been  shown 
to  have  both  commercial  and  military  applications  for 
obtaining  significant  inclination  changes  for  aerodynamically 
assisted  orbital  transfers.  The  improved  efficiency  of  the 
aerobang  technique  over  the  aerocruise  for  subcircular  and 
supercircular  speeds  at  two  discrete  heating  rate  values  given 
in  Reference  3  was  shown  here  to  extend  to  the  entire  range  of 
heating  rates  studied,  not  just  the  high  and  low  heating  rate 
cases . 

Thrust  vector  control  can  be  effective  in  increasing 
inclination  changes  for  both  short  and  long  duration  aerobang 
maneuvers.  This  was  demonstrated  for  both  supercircular  and 
subcircular  velocities.  The  use  of  thrust  vector  control  was 
most  effective  when  applied  at  a  constant  angle  of  attack. 

The  effects  of  increasing  mass  fraction  of  fuel  consumed 
from  two  to  twenty  percent  had  a  significant  effect  on  the 
aerobang  maneuver,  beyond  the  obvious  extension  of  the  time 
used  to  perform  the  maneuver.  For  the  subcircular  maneuvers, 
the  constant  heating  rate  method  of  controlling  the  aerobang 
flight  resulted  in  excellent  inclination  changes  but  with 
extremely  high  (70°)  angles  of  attack.  By  controlling  the 
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angles  of  attack,  the  heating  rates  increase  by  nearly  a 
factor  of  two,  placing  new  requirements  on  any  spacecraft 
design . 

The  supercircular  aerobang  maneuvers  were  also  simulated 
for  both  constant  heating  rate  and  constant  angle  of  attack 
cases.  The  constant  heating  rate  case  proved  to  be 
ineffective  in  generating  even  one  degree  of  inclination 
change  because  of  the  loss  of  angle  of  attack  and,  therefore, 
lift.  The  constant  angle  of  attack  runs,  however,  proved  to 
be  almost  as  effective  as  for  the  subcircular  cases.  This 
case  also  provided  two  inherent  advantages  over  the 
subcircular  case,  an  increasing  altitude  which  translates  to 
a  decreasing  heating  rate  value  and  the  likelihood  of  an 
easier  reorbit  maneuver. 

B.  RECOMMENDATIONS 

There  remain  several  areas  in  the  study  of  aerodynamically 
assisted  orbital  plane  change  maneuvers  that  have  yet  to  be 
explored.  Among  these  are  the  effects  of  combining  any  two  or 
all  of  the  three  techniques  (aeroglide,  aerocruise,  and 
aerobang)  to  maximize  the  fuel  efficiency  of  a  given 
inclination  change  requirement. 

This  thesis,  as  well  as  most  other  work  in  this  area, 
deals  only  with  the  atmospheric  passage  portion  of  the  orbit 
transfer.  By  extending  the  simulations  to  cover  an  entire 
orbit  transfer,  including  both  entering  and  exiting  the 
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atmosphere,  a  better  measure  of  the  effectiveness  of  these 
maneuvers  can  be  made. 


This  is  the  first  work  in  which  the  constant  angle  of 
attack  method  of  controlling  the  aerobang  maneuver  was  used. 
Further  studies  should  concentrate  on  whether  or  not  it  can 
indeed  be  a  viable  means  of  control,  from  both  a  flight 
control  system  and  thermal  protection  point  of  view. 

In  order  to  verify  the  applicability  of  aerodynamical ly 
assisted  maneuvers,  the  study  of  thermal  protection  system 
requirements  and  techniques  must  be  included.  Superior 
atmospheric  and  heating  models  can  also  be  used  to  refine  the 
requirements  for  such  systems. 

Finally,  a  means  of  minimizing  the  time  required  for  a 
maneuver  to  meet  a  given  inclination  change  requirement  could 
be  established.  This  would  allow  the  mission  goal  to  be  met 
while  minimizing  the  integrated  heat  load  on  the  spacecraft. 
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PROGRAM  ORBIT 
IMPLICIT  REAL*8(A-H,L-Z) 

DIMENSION  X(6) ,XDOT(6) ,C(6) 

PROGRAMMER:  THOMAS  P.  SPRIESTERBACH 
MODIFIED  BY:  RICHARD  E.  JOHNSON 
DATE  :  15  FEB  93 

SYSTEM  :  IRIS  INDIGO/FORTRAN 


VARIABLE  LIST 


I 

INDEX 

J 

KOUNT 

DUM 

OLDCO 

CHANGE 

NE 

T 

TF 

H 

TP  I 

X(l) 

X(2) 

X(3) 

X(4) 

X(5) 

X(6) 

XDOT( ) 

C(l) 

C(2) 

C(3) 

C(4) 

C(5) 

C(6) 

BETA 

GO 

G 

HO 

RHOO 

MU 


COUNTING  INDEX 

COUNTING  INDEX  FOR  RUNGE-KUTTA  ROUTINE 
COUNTING  INDEX 

COUNTING  INDEX  FOR  OUTPUT  DETERMINATION 
DUMMY  VARIABLE  FOR  CONVENIENCE 
CONTROL  VARIABLE  FOR  CONVERGENCE  ON  ANGLE  OF 
ATTACK  [AOA] 

INCREMENTAL  CHANGE  IN  AOA  USING  NEWTONS  METHOD 


NUMBER  OF  EQUATIONS  OF  MOTION 
TIME  (SEC) 

FINAL  TIME  (SEC) 

INCREMENT  OF  TIME  (RUNGE-KUTTA  ROUTINE]  (SEC) 

TIME  PRINT  INTERVAL  (SEC) 

RADIUS  (METERS) 

THETA  (RADIANS) 

PHI  (RADIANS) 

VELOCITY  (M/SEC) 

FLIGHT  PATH  ANGLE  GAMMA  (RADIANS) 

PSI  (RADIANS) 


DERIVATIVES  OF  THE  ABOVE  SIX  VARIABLES 
ZEROTH  ORDER  COEFFICIENT  FOR  CL  EQUATION 
FIRST  ORDER  COEFFICIENT  FOR  CL  EQUATION 
SECOND  ORDER  COEFFICIENT  FOR  CL  EQUATION 
ZEROTH  ORDER  COEFFICIENT  FOR  CD  EQUATION 
FIRST  ORDER  COEFFICIENT  FOR  CD  EQUATION 
SECOND  ORDER  COEFFICIENT  FOR  CD  EQUATION 
DENSITY  MODEL  EXPONENT  (METERS) 

GRAVITY  AT  EARTH  SURFACE  (  9 . 806M''2/S ) 

LOCAL  GRAVIATIONAL  ATTRACTION  (M'-Z/S) 

REFERENCE  ALTITUDE  FOR  DENSITY  MODEL  (METERS) 

REFERENCE  DENSITY  FOR  DENSITY  MODEL  (KG/M^3) 

GRAVITATIONAL  CONSTANT 
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c 

ALPHA 

ANGLE  OF  ATTACK 

(RADIANS) 

c 

ALFA 

ANGLE  OF  ATTACK 

(DEG) 

c 

VEC 

THRUST  VECTORING  ANGLE 

(DEG) 

c 

EPS 

THRUST  VECTORING  ANGLE 

(RADIANS) 

c 

MASS 

MASS  OF  SPACECRAFT 

(KG) 

c 

MASSO 

INITIAL  MASS  OF  SPACECRAFT 

(KG) 

c 

SPI 

SPECIFIC  IMPULSE  OF  POWER  PLANT 

(  SEC) 

c 

R 

RADIUS 

(M) 

c 

V 

VELOCITY 

(M/S) 

c 

RADIUS 

RADIUS 

(KM) 

c 

VELOCITY 

VELOCITY 

(KM/S) 

c 

GAMMA 

FLIGHT  PATH  ANGLE 

(DEG) 

c 

ANGM 

ANGULAR  MOMENTUM 

(M"2/S) 

c 

SIGMA 

ANGLE  OF  BANK 

(DEG) 

c 

OINC 

INCLINATION 

(DEG) 

c 

DELTA 

INTERMEDIATE  ANGLE 

(DEG) 

c 

OMEGA 

LONGITUDE  OF  THE  ASCENDING  NODE 

(DEG) 

c 

AS 

TANGENTIAL  ACCELERATION 

(M/S-2) 

c 

AR 

RADIAL  ACCELERATION 

(M/S-2) 

c 

AW 

BINORMAL  ACCELERATION 

(M/S'2) 

c 

THR 

THRUST 

(N) 

c 

AOB 

ANGLE  OF  BANK 

(RADIANS) 

c 

RHO 

DENSITY 

(KG/M'3 ) 

c 

S 

REFERENCE  AREA  OF  SPACECRAFT 

(M-2) 

c 

N 

DENSITY  COEFFICIENT  FOR  HEATING  RATE 

EQUATION 

c 

M 

VELOCITY  COEFFICIENT  FOR  HEATING  RATE  EQUATION 

c 

HEAT 

STAGNATION  POINT  HEATING  RATE 

(WATTS/M'-2  ) 

c 

QS 

DYNAMIC  PRESSURE  TIMES  THE  REFERENCE 

AREA  (N*M"2) 

c 

CD 

COEFFICIENT  OF  DRAG 

c 

CL 

COEFFICIENT  OF  LIGT 

c 

LIFT 

LIFT  ON  SPACECRAFT 

(N) 

c 

DRAG 

DRAG  ON  SPACECRAFT 

(N) 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccrccccccccccccccc 

c 

C  DEFINED  FUNCTION  FOR  ACCEL  DUE  TO  GRAVITY 

C 

G(R)=MU/(R*R) 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 

C  OPEN  TWO  INPUT  FILES . DATA.DAT  -  ORBITAL  PARMS 

C  AERO.DAT  -  VEHICLE  PARMS 

C  ONE  OUTPUT  FILE . OUT. DAT  -  TABLE  OF  VALUES 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 

OPEN(10,FILE='data.dat' ,  STATUS='OLD' ) 

OPEN( 12,FILE='aero.dat ' , STATUS='OLD' ) 
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OPEN( 13,FILE='out.dat' ,  STATUS='NEW' ) 


READ( 10,1)(X(I),I=1,6),T,TF,H,TPI 

1  FORMAT{/,/,20X,10(/,20X,Dl3.7)) 

READ( 12,2)(C(I),I=1,6) , AOPT, AOB , N , M, S , MASSO , FM , SP I , THR , EPS 

2  FORMAT(/, /,/,/, 16(/,20X,D13. 7) ) 

C 

C  NUMBER  OF  EQUATIONS  TO  INTERGRATE 

C  X( 1 )... .RADIUS 

C  X(2) _ THETA  (SPHERICAL  COORD.  PARAMETER) 

C  X(3) - PHI  (SPHERICAL  COORD.  PARAMETER) 

C  X(4) ... .VELOCITY 

C  X(5) ... .FLIGHT  PATH  ANGLE  GAMMA 

C  X(6)....PSI 

C 

C  INITIALIZATION  OF  CONSTANTS 

C 

NE  =  6 

PI=3.14169265359D+0 
MU=3.986012D+14 
G0=9.806D+0 
MASS=MASSO 
VEC=EPS*180./PI 

MAIN  PROGRAM 

INDEX=0 
KOUNT=l 

INITIAL  CALL  TO  GET  ACCELERATIONS  AND  CONTROL  VARIABLES 
FOR  FIRST  OUTPUT 

CALL  CNTRL ( DRAG , L I FT , THR , ALPHA , AOB , X , MASS , S , C , N , M , RHO , EPS ) 
CALL  ACEL ( AS , AR , AW , DRAG , L I FT , THR , EP S , ALPHA , AOB , MAS S ) 

CALL  HDR(X( 1) ,X(4) ,MASS,THR,VEC,X(5) ) 

CALL  WRT(X,T,THR,MASS, ALPHA, AOB, RHO, N,M) 


THIS  IS  THE  MAIN  BLOCK  OF  THE  PROGRAM  CALLING  FIRST  THE 
CONTROL  SUBROUTINE  TO  DETERMINE  THE  OUTPUT  VARIABLES.  NEXT, 
THE  ACCELERATION  ROUTINE  CALCULATES  THE  RELATIVE 
ACCELERATION  TO  THE  VEHICLE,  THEN  THE  DIFFERENTIAL  EQUATIONS 
OF  MOTION  ARE  DEFINED  AND  A  FOURTH  ORDER  RUNGE-KUTTA  ROUTINE 
IS  USED  TO  INTEGRATE  THEM.  THIS  BLOCK  IS  COMPUTED  AT  EACH 
TIME  INCREMENT  H. 


•  •  •  • 


•  •  • 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


100 


CALL  CNTRL ( DRAG , L I FT , THR , ALPHA , AOB , X , MAS S , S , C , N , M , RHO , EPS ) 
CALL  ACEL ( AS , AR , AW , DRAG , L I FT , THR , EPS , ALPHA , AOB , MAS S ) 

CALL  ORB ( X , XDOT , AS , AR , AW , G ( X ( 1 ) ) ) 

CALL  RK4 (T,X, XDOT, NE,H, INDEX) 

IF (INDEX  .NE.  0)  GO  TO  100 


Ar 


C  DETERMINE  IF  TIME  TO  PRINT  TO  OUTPUT  USING  TIME  PRINT 

C  INTERVAL 


200  IF(KOUNT  .LT.  IDNINT ( TPI /H ) )  GO  TO  300 

CALL  WRT ( X , T , THR , MAS S , ALPHA , AOB , RHO , N , M ) 
C  PRINT  STATUS  TO  TERMINAL 


WRITE ( 6, 7 )T, ALPHA, AOB 
7  FORMAT( ' ,F8. 1 , IX, F6 . 3, 1X,F6 . 3) 

KOUNT=0 

300  KOUNT=KOUNT+l 
C 

C  CHANGE  MASS  OF  S/C 

C 

MASS=MASS-THR*H/( SPI*G0 ) 

C 

C  TWO  TESTS  EITHER  CAN  STOP  THE  PROGRAM.  TEST  FOR  MASS  LESS 

C  THAN  FINAL  MASS  AND  ALSO  IF  TIME  IS  GREATER  THAN  FINAL  TIME 

C 

IF(MASS  .GT.  FM  .AND.  T  . LE .  TF)  GO  TO  100 
C 

C  PRINT  OUTPUT  AGAIN  IF  MASS  SWITCH  ENDED  PROGRAM 

C 

IF(MASS  .LE.  FM)  CALL  WRT ( X , T , THR, MASS , ALPHA , AOB , RHO , N , M ) 

C 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  THIS  SUBROUTINE  CALCULATES  THE  CONTROL  OF  THE  FREE  VARIABLES 

C  DURING  THE  ATMOSHPHERIC  FLIGHT. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE  CNTRL ( DRAG , LIFT , THR , ALPHA , AOB , X , MASS , S , 
*C,N,M,RHO,EPS) 

IMPLICIT  REAL*8(A-H,L-Z) 

DIMENSION  X(6),C(6) 

C 

MU=3.986012D+14 
GAMMA=X ( 5 ) 


•  • 


•  • 


•  • 
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nno  nnnn  nnnnnn  ooonnn  nnnonono  nnnno 


R=X(1) 

V=X(4) 

COEFFICIENTS  FOR  AN  EXPONENTIAL  DENSITY  MODEL 
REFERENCE  J  MEASE  1976  US  STANDARD  ATMOSPHERE 
RANGE  50KM  TO  nOKM 

BETA=1 .41D-4 

RHO0=3.0968D-4 

R0=6,435D+6 

G=MU/(R*R) 

RHO=RHOO*EXP(-BETA*(R-RO) ' 

QS= . 5D+0*RHO*V*V*S 


FOR  AEROBANG  THE  AOB  IS  SET  AND  THE  ANGLE  OF  ATTACK  IS 
CONTROLLED  TO  FLY  AT  A  CONSTANT  HEATING  RATE;  THE  ALTITUDE 
AND  VELOCITY  ARE  ALLOWED  TO  FLOAT. 


GUESS  ALPHA  INITIALLY  FOR  NEWTON  APPROXIMATION  ROUTINE 

ALPHA= ( C ( 5 ) *QS  +  SQRT ( C ( 5 ) *  *2  +  4 . *QS*C ( 6 ) * ( THR*COS ( EPS ) 

*  -QS*C(4) ) ) )/(2.*C(6)*QS) 

USE  NEWTON  APPROXIMATION  METHOD  TO  CONVERGE  ON  ALPHA. 

USE  A  WHILE  STRUCTURE  FOR  THE  CONVERGENCE  OF  ALPHA. 

SET  OLDCO  AND  CHANGE  EQUAL  TO  THE  VALUES  BELOW  TO  ENSURE 
THAT  THE  PROGRAM  SWITCHES  ARE  NOT  INITIALLY  TRIPPED. 

OLDCO=l .D+10 
CHANGE=1 . 

OR,  ALLOW  ALPHA  TO  BE  A  CONSTANT  AND  OBSERVE  THE  BEHAVIOR  OF 
THE  HEATING  RATE  AND  INCLINATION. 

ALPHA  =  0.69813 

DO  WHILE  (ABS(CHANGE)  .GT.  l.D-3) 

CL=C ( I ) +C ( 2 ) * ALPHA+C ( 3 ) *ALPHA*  *2 
CD=C ( 4 ) +C ( 5 ) *ALPHA+C ( 6 ) *ALPHA**2 

CALCULATE  THE  LIFT  AND  DRAG  FORCES  FROM  THE  COEFFICIENTS 
CALCULATED  ABOVE. 

LIFT=CL*QS 

DRAG=CD*QS 

THESE  THREE  EQUATIONS  REPRESENT  THE  FUNCTION  AND  ITS 
DERIVATIVE  WITH  RESPECT  TO  ALPHA  FOR  NEWTON'S  METHOD  OF 
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C  APPROXIMATION 

C 

COEFF= ( THR*COS ( ALPHA+EPS ) -DRAG ) -MASS*  S IN ( GAMMA ) 

** (G+BETA* (N/M) *V*V) 

COEFP=-THR*SIN( ALPHA+EPS) - (C( 5 )+2*C( 6 ) *ABS( ALPHA) ) *QS 
C 

CHANGE=COEFF/COEFP 

ALPHA=ALPHA-CHANGE 

IF  ( ABS(OLDCO)-ABS( CHANGE ) .LE.O. )  STOP 
OLDCO=CHANGE 
END  DO 
RETURN 
END 
C 
C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  THIS  SUBROUTINE  DEFINES  THE  TANGENTIAL,  NORMAL,  AND  BINORMAL 

C  ACCELERATIONS  WITH  RESPECT  TO  THE  VEHICLE  AND  THE  ORBITAL 

C  PLANE . 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  ACEL( AS, AR, AW, DRAG, LIFT, THR, EPS , ALPHA , AOB , MASS ) 
IMPLICIT  REAL*8(A-H,L-Z) 

C 

C  COMPUTE  THE  ACCELERATIONS  ON  THE  FLIGHT  VEHICLE  USING 

C  THE  THREE  ACCELERATION  EQUATIONS,  WHICH  INCLUDE  THE 

C  THRUST  VECTORING  ANGLES. 

C 

AS=(THR*COS( ALPHA+EPS) -DRAG) /MASS 
AW= ( LI FT+THR* SIN (ALPHA+EPS) ) * S IN ( AOB ) /MASS 
AR=(LIFT+THR* SIN (ALPHA+EPS) ) *COS ( AOB ) /MASS 
C 

RETURN 

END 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C  THIS  SUBROUTINE  IS  THE  COLLECTION  OF  DIFFERENTIAL  EQUATIONS 

C  WHICH  DESCRIBE  THE  MOTION  OF  THE  SPACECRAFT. 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE  ORB ( X , XDOT , AS , AR , AW, G ) 

IMPLICIT  REAL*8(A-H,L-Z) 

DIMENSION  X(6),XDOT(6) 

C 

XDOT(l)=X(4)*SIN(X(5) ) 


XD0T(2 )=X(4 )*COS(X( 5) ) *C0S(X(6) )/(X( 1 ) *COS(X( 3 ) ) ) 
XDOT(3)=X(4)*COS(X(5) ) *SIN(X( 6 ) ) /X( 1 ) 

XDOT(4)=AS-G*SIN(X(5) ) 

XDOT( 5 )=( AR-G*COS(X( 5 ) ) +X( 4 ) *X( 4 ) *COS(X( 5 ) ) /X ( 1 ) ) /X ( 4 ) 
DUM=TAN(X(3) )/X(l) 

XDOT(6)=AW/(COS(X(5) )*X(4) ) -X ( 4 ) *COS ( X ( 5 ) )*C0S(X(6) )*DUM 
C 

RETURN 

END 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  THIS  SUBROUTINE  TAKES  CARE  OF  ALL  HARDCOPY  OUTPUT 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE  WRT ( X , T , THR , MASS , ALPHA , AOB , RHO , N , M ) 

IMPLICIT  REAL*8(A-H,L-Z) 

DIMENSION  X(6),XDOT(6) 

PI=3 . 14159265359D+0 

RADIUS=X( 1 )/1000. 

VELOCITY=X(4)/1000, 

GAMMA=X(5)*180./PI 
ALFA=ALPHA*180, /PI 
HEAT=RHO**N*X( 4 ) **M*9 . 652D-5 
SIGMA=AOB*180,/PI 
OINC=ACOSD(COS(X(3) )*COS(X(6) ) ) 

DELTA=0.0 

IF  (OINC  .NE.  0.)  DELTA=ASIND(TAN(X(3) )/TAND(OINC) ) 
OMEGA=X(2 ) *180 . /PI-DELTA 

WRITE ( 13, 1 )T, RADIUS, VELOCITY, MASS, GAMMA, OINC, THR, HEAT, 
*ALFA, SIGMA 

1  FORMAT { IX, F6. 1, 1X,F8.3, 1X,F6.4,2X,F7 ,2, IX, F5 . 3, IX, 

*  F6 . 3, 1X,F7 . 1 , 1X,E9 . 3, 1X,F7 . 3, 1X,F7. 3 ) 

C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  THIS  SUBROUTINE  ATTACH  A  HEADER  TC  THE  OUTPUTS  TO  KEEP  TRACK 

C  OF  THE  DIFFERENT  CHANGES 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE  HDR(R, V, MASS, THR, VEC, GAMMA) 

IMPLICIT  REAL*8(A-H,L-Z) 


noon 


DO  1  1=13,14 

IF  ( 1 .EQ. 14 )GO  TO  1 

WR1TE( 1 , * ) 

WR1TE(1,*)'  SELECTED  INITIAL  INPUT  DATA:' 

WRITE ( 1, * ) 

WRITE ( I,2)R,V,MASS,THR,VEC,GAMMA 
WRITE( I, *) 

WRITE (I,*) 

CONTINUE 

FORMAT ( '  INITIAL  RADIUS  (METERS)  ',F10.2,/,'  INITIAL 

♦VELOCITY  (METERS/SEC)  ' , F7 . 2 , / , '  INITIAL  MASS  (KG) 

♦  ',F7.2,/,'  INITIAL  THRUST  (NEWTONS) 

*,F9.2,/,'  THRUST  VECTOR  ANGLE  (DEG) 

*,F6.2,/,'  INITIAL  FLIGH":  PATH  ANGLE  (RADS)  '  ,  F5 . 2  ) 


WRITE(13,*)'  TIME  RADIUS  VELOCITY  MASS  GAMMA  INCLI 


*  THRUST  QDOT 
WRITE(13,*)'  (SEC) 

*  (N)  (J/M2S) 


ALPHA  AOB' 

(KM)  (KM/SEC)  (KG)  (DEG)  (DEG) 
(DEG)  (DEG)' 


RETURN 

END 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

THIS  SUBROUTINE  IS  A  FOURTH  ORDER  RUNGE-KUTTA  ROUTINE  TAKEN 
FROM  PROF.  ROSS'S  SUBROUTINE. 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE  RK4 ( T , X , XDOT , NE , H , INDEX) 

IMPLICIT  REAL*8(A-H,L-Z) 

INTEGER  INDEX, I 

DIMENSION  X( 6) ,XDOT( 6) , SAVED(6) , SAVEX( 6) 

C 

INDEX=INDEX+1 

GO  TO  (  1,2, 3,4) , INDEX 

1  DO  10  1=1, NE 
SAVEX( I )=X( I ) 

SAVED ( I)=XDOT( I) 

10  X(I)=SAVEX(I)+.5DO*H*XDOT(I) 

T=T+.5D0*H 

RETURN 

C 

2  DO  20  1=1, NE 

SAVED( I )=SAVED( I)+2 .DO*XDOT( I) 

20  X(I)=SAVEX(I)+.5DO*H*XDOT(I) 

RETURN 


DO  30  1=1, NE 
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u 


30 


SAVED  (  I )  =SAVED(  I )  -*-2  .  DO  *XDOT  (  I ) 

X( I ) =SAVEX( I ) +H*XDOT( I ) 

T=T+.5D0*H  • 

RETURN 

DO  40  I=1,NE 
40  X(  I  )=SAVEX(  I  )  -fH/6  .  DO*  (SAVED (  I )  +XDOT(  I )  ) 

INDEX=0 

RETURN  • 

END 
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•  •  • 


•  •  •  •  • 


SAMPLE  AERO.DAT  INPUT  FILE 


This  file  gives  the  aerodynamic  data  used  for  the  MRRV. 
Cl-6, AOPT,N,M,S,MASS,FM, SP I, THRUST 


COEFFICENT  Cl 
COEFFICENT  C2 
COEFFICENT  C3 
COEFFICENT  C4 
COEFFICENT  C5 
COEFFICENT  C6 
AOA  FOR  CL/CD  MAX 
AOB  FOR  BANG  MANU 
DENSITY  EXPONENT 
VELOCITY  EXPONENT 
REFERNCE  AREA 
INITIAL  MASS  (KG) 
FINAL  MASS  (KG) 
SPECIFIC  IMPULSE  (S) 
THRUST  (N) 

THRUST  ANGLE  (RAD) 


l.OOOOOOD-02 
2.860000D-01 
1 . 313000D+00 
4.700000D-02 
4.470000D-01 
2 . 040000D-t-00 
1 .823870D-01 
1 .570796D-t-00 
0.500000D+00 
3. 150000D+00 
1 . 169800D+01 
4.898000D+03 
4.800000D+03 
2.950000D+02 
1 .467900D+04 
0.261800D+00 


SAMPLE  DATA. OAT  INPUT  FILE 


This  file  gives  the  orbital  and  time  printing  information  for 
the  simulation  program. 


INPUT  ORBITAL  PARAMETERS 


RADIUS  METERS 
THETA  RADIANS 
PHI  RADIANS 
VELOCITY  M/S 
FLIGHT  PATH  ANGLE 
PSI  RADIANS 
BEGIN  TIME  SEC 
END  TIME  SEC 
TIME  INTERVAL  SEC 
PRINT  TIME  INTERVAL 


6.4450000D-t-06 
O.OOOOOOOD+00 
0. OOOOOOOD+03 
7 . 7100000D+03 
O.OOOOOOOD+00 
O.OOOOOOOD+00 
O.OOOOOOOD+00 
4.0000000D-t-03 
5.0000000D-03 
1 . OOOOOOOD+00 
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•  •  • 


•  • 


SAMPLE  OUT. OAT  OUTPUT  FILE 


SELECTED  INITIAL  INPUT  DATA: 


INITIAL  RADIUS  (METERS) 

INITIAL  VELOCITY  (METERS/SEC) 

INITIAL  MASS  (KG) 

INITIAL  THRUST  (NEWTONS) 

THRUST  VECTOR  ANGLE  (DEG) 

INITIAL  FLIGHT  PATH  ANGLE  (RADS) 


6445000.00 

7710.00 

4898.00 

14679.00 

15.00 

0.00 


TIME 

RADIUS  ' 

VELOCITY 

MASS 

GAMMA 

INCLI 

THRUST 

QDOT 

ALPHA 

AOB 

SEC) 

(KM) 

(KM/SEC) 

(KG) 

(DEG) 

(DEG) 

(N) 

( J/M2S ) 

(DEG) 

(DEG) 

0.0 

6445.000 

7.7100 

4898 . 00 

0.000 

0.000 

14679.0 

0. 147E  +  07 

30.615 

90.000 

1.0 

6445.000 

7.7100 

4892.95 

-  .003 

0.037 

14679.0 

0  .  147E  +  07 

30.937 

90.000 

2.0 

6444 . 999 

7 .7099 

4887 . 88 

-  .  006 

0.074 

14679.0 

0.  147E  +  07 

31.254 

90.000 

3.0 

6444.998 

7.7097 

4882.80 

-.008 

0.112 

14679.0 

0. 147E  +  07 

31.567 

90.000 

4.0 

6444 . 997 

7.7095 

4877.73 

-.011 

0.150 

14679.0 

0. 147E  +  07 

31.876 

90.000 

5.0 

6444.995 

7.7092 

4872 . 65 

-.014 

0.189 

14679.0 

0. 147E  +  07 

32 . 181 

90.000 

6.0 

6444.993 

7.7088 

4867 . 58 

-  .017 

0.228 

14679.0 

0. 147E  +  07 

32.482 

90.000 

7.0 

6444.991 

7.7084 

4862 . 50 

-.019 

0.268 

14679.0 

0. 147E  +  07 

32.779 

90.000 

8.0 

6444.988 

7.7079 

4857.43 

-.022 

0.309 

14679.0 

0. 147E+07 

33.074 

90.000 

9.0 

6444.985 

7 . 7074 

4852.36 

-.025 

0.349 

14679.0 

0. 147E  +  07 

33.365 

90 . 000 

10.0 

6444.981 

7.7068 

4847.28 

-  .028 

0.391 

14679.0 

0. 147E  +  07 

33.653 

90.000 

11.0 

6444.977 

7.7061 

4842.21 

-.031 

0.433 

14679.0 

0. 147E+07 

33.938 

90.000 

12.0 

6444.973 

7.7053 

4837 .13 

-.034 

0.475 

14679.0 

0. 147E+07 

34.221 

90.000 

13.0 

6444.968 

7.7045 

4832.06 

-.036 

0.518 

14679.0 

0. 147E+07 

34.501 

90.000 

14.0 

6444.963 

7 . 7037 

4826.98 

-  .039 

0.561 

14679.0 

0. 147E  +  07 

34.778 

90.000 

15.0 

6444.958 

7 .7027 

4821.91 

-  .042 

0.605 

14679.0 

0. 147E  +  07 

35.054 

90.000 

16.0 

6444.952 

7.7017 

4816 . 84 

-  .045 

0.649 

14679.0 

0. 147E  +  07 

35.327 

90.000 

17.0 

6444.946 

7 .7006 

4811.76 

-.048 

0.694 

14679.0 

0 . 147E  +  07 

35.598 

90.000 

18.0 

6444.939 

7.6995 

4806.69 

-.051 

0.739 

14679.0 

0. 147E  +  07 

35.868 

90.000 

19.0 

6444.932 

7.6983 

4801.61 

-.054 

0.785 

14679.0 

0 . 147E  +  07 

36.135 

90.000 

19.3 

6444.930 

7.6979 

4799.99 

-.055 

0.800 

14679.0 

0. 147E  +  07 

36.219 

90.000 
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